Dynamics of the BCS-BEC crossover in a degenerate Fermi gas 
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We study the short-time dynamics of a degenerate Fermi gas positioned near a Feshbach resonance 
following an abrupt jump in the atomic interaction resulting from a change of external magnetic 
field. We investigate the dynamics of the condensate order parameter and pair wavefunction for a 
range of field strengths. When the abrupt jump is sufficient to span the BCS to BEC crossover, 
we show that the rigidity of the momentum distribution precludes any atom-molecule oscillations 
in the entrance channel dominated resonances observed in the 40 K and 6 Li. Focusing on material 
parameters tailored to the 40 K Feshbach resonance system at 202.1 gauss, we comment on the 
integrity of the fast sweet projection technique as a vehicle to explore the condensed phase in the 
crossover region. 
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Ultracold alkali atomic gases provide a valuable arena 
in which to explore molecular Bose-Einstein condensa- 
tion (BEC) 1] and fermionic pair condensation 0, 0- 
Further, the facility to control interparticle interactions 
via a magnetically-tuned Feshbach resonance (FR) pro- 
vide a unique opportunity to investigate the BCS-BEC 
crossover and the dynamics of condensate formation. As 
well as the adiabatic association of molecules pj, both 
fast sweep 'projections' of fermionic pair condensates 
onto the molecular BEC 0, Q , and atom-molecule Ram- 
sey fringes |^ have been reported in the recent literature. 
Lately, motivated by earlier work on the mean-field BCS 
system [f|, it was shown in separate works 0, Q that 
the mean-field equations of motion of a Bose-Fermi (BF) 
model, commonly used to describe the FR system, are 
characterised by an integrable nonlinear dynamics. From 
these works, three striking predictions emerged: Firstly, 
when perturbed by an abrupt change in the strength 
of the pair interaction, the condensate order parameter 
exhibits substantial oscillations which range in magni- 
tude between some initial state value, At, and that ex- 
pected for the equilibrium final state configuration A eq . 
Secondly, in the absence of energy relaxation processes, 
these oscillations remain undamped suggesting the po- 
tential to observe coherent atom-molecule oscillations in 
the FR system. Thirdly, it is proposed that a spectral 
'hole-burning' phenomena in the atomic momentum dis- 
tribution, at a frequency associated with one half of the 
molecular binding energy, provides a signature of such 
atom-molecule oscillations 0. 

In the following, we will argue that this behaviour rests 
on an auxiliary constraint that, in the present system, 
seems hard to justify. Drawing on the results of a nu- 
merical analysis of the unconstrained dynamics, we will 
show that, in the absence of relaxational processes, the 
oscillations of the order parameter are damped substan- 
tially, even at the level of the mean-field. When the 
abrupt change of the interaction is sufficient to span the 
BCS-BEC crossover, for both the single channel and BF 
systems, the magnitude of oscillations remain small at- 
tenuating to some value Ap which lies close to the initial 



state, and much smaller than the expected equilibrium 
state value, A oq , while the distribution remains essen- 
tially "frozen" to the initial state. The hole-burning os- 
cillations of the form predicted in Ref. Q do not appear 
in either regime. From these results, we are able to con- 
clude that, although the observation of BCS-BEC like 
atom-molecule oscillations in alkali Fermi gases seems in- 
feasible, the rigidity of the initial state distributions val- 
idate fast sweep techniques as a reliable experimental 
tool to explore the crossover region. 

Formally, a detailed microscopic theory of FR phenom- 
ena demands consideration of all matrix elements con- 
necting different spin states participating in condensate 
formation. In practice, an accurate description of the res- 
onance can be obtained either by using a single-channel 
magnetic field dependent effective interaction between 
atoms in the entrance channel or, more widely, in the two- 
most relevant channels. For low relative momenta, rele- 
vant to cold atom physics, the full form of complex atomic 
potentials are not resolved. Indeed, separable potentials 
with microscopic parameters drawn from experiment and 
exact multichannel calculations can be used to recover all 
low-energy binary scattering observables 0,0. Since, 
very often, only one bound state of the closed channel po- 
tential is relevant, it has been traditionally replaced by a 
fictitious Bose-particle and FR phenomena captured by 
an effective two-channel BF Hamiltonian [i^ll3|. 

Now, applied to the entrance channel dominated reso- 
nances observed in the 40 K and 6 Li system, formal calcu- 
lations [Tol[lll| support a picture in which the BCS-BEC 
crossover is mediated by only a very small admixture of 
closed channel states — e.g., in the K FR at 202.1 
gauss, the admixture of the closed channel is ca. 8% or 
less of the total ^]J. Rather, the weakly bound molec- 
ular state appears at a detuning Eq which lies far from 
the value where the resonance state of the closed channel 
crosses the dissociation threshold. Since the two-body 
observables drawn from the exact numerical solution of 
the Shrodinger equation within finite-range single and 
two-channel models do not differ over a wide range of 
fields j FR phenomena in 40 K can be equally- well de- 
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scribed by a single-channel theory, 
H = ^2e k al s a Us +Y^ ^'4+qT a -k4 a -k'.l. k'+q'r' (1) 

ks kk'q 

involving Fermi operators a ks and a ks . In the following, 
to account for the entire region of crossover (and not only 
the universal regime) we take as matrix elements Vkk' = 
Vb(-S)xk(ffbg)Xk'(c r bg) with Xk^bg) = exp[-(kcr bg ) 2 /2] 
and the parameters Vq (B) and Obg chosen to recover the 
correct magnetic field dependence of the scattering length 
and the highest vibrational bound state 0] . 

Applied to the fields $ k = X} s ( a iL aks )/2 and K k = 
(a_kj,dk,|), the Heisenberg equations of motion i{A) = 
([A, H]) translate to the relations 

ik u = 2e k K k - A k (2$ k - 1) 

= A k K k -A k K k , (2) 

where A k = Vo(B)xu.(crbg) X)k' Xk'(oiag)ftk' denotes the 
complex order parameter. Similarly for a BF theory, 
defining g k = .goXfe(c) as the coupling of the entrance 
channel states to the bosonic field 6 k associated with the 
Feshbach resonance level configuration of the closed chan- 
nel, the equations of motion acquire the same form as @ 
with A k = 3oXk(cr)ko + VbgXk(ffbg) J2w Xk'Obg)Kk' and 
fro = (frk=o) obeying the supplementary equation ibo = 
Eo(B)bo + 5o X)k Xk(c)K k . As with the single-channel 
theory, the five parameters which characterise the reso- 
nance, the background potential strength Vb g and range 
Cbg (which define the entrance channel scattering length 
and its highest vibrational bound state), the interchan- 
nel coupling go and its range a, and the detuning Eq(B) 
(which specifies the position and width of the resonance), 
are determined from experiment and exact multichannel 
calculations [TT| . 

Before turning to the numerical analysis, it is instruc- 
tive to contrast the present approach to that adopted 
in Refs. [5J, |fj [jj • Given an initial condition, the Heisen- 
berg equations of motion (J2J) present a deterministic time 
evolution of the density distributions. Indeed, the con- 
servation of total atomic density n, implicit in the dy- 
namics |PJ, provides a check on the integrity of the nu- 
merical integration described below. By contrast, the 
integrability of the equations of motion as described by 
Refs. |5|,|6|jl3 relies on an additional constraint involving 
the density and an auxiliary parameter playing the role of 
a "chemical potential" . The need for this constraint orig- 
inates from the fact that, as with the equilibrium steady 
state solution, the analytical solution of the dynamical 
equations is determined only up to an arbitrary momen- 
tum dependent sign factor which is fixed by the auxiliary 
constraint pj , instead of naturally following from the ini- 
tial conditions. A similar phenomenology describes the 
effect of a classical laser source on a semiconductor elec- 
tron/hole system where the laser frequency "imprints" 
a chemical potential onto the system While equili- 
bration processes are sufficiently small, the electrons and 



holes assume a non-equilibrium distribution around the 
externally imposed chemical potential resulting in a phe- 
nomenon of "spectral-hole burning" in which the density 
distribution is depleted at the laser frequency. In the 
atomic gas system, where the degrees of freedom remain 
internal, it is difficult to see how such a choice is mo- 
tivated or justified. Crucially, we will see that the un- 
constrained dynamics associated with the equations of 
motion J5J) lead to physical behaviour very different from 
that obtained from the constrained dynamics 0, IE • 

With this background, let us now turn to the results 
of the numerical investigation of the dynamics @ per- 
formed using an adaptive step Runge-Kutta algorithm 
applied on a grid in momentum space chosen fine enough 
to ensure convergence at each time step. Although we 
find that the qualitative behaviour of BCS and BF dy- 
namics is generic, we focus specifically on potentials char- 
acteristic for the 40 K resonance at Bq = 202 .1 gauss with 
a density of n = 1.5 x 10 13 cm" 3 (i.e. T F = 0.35 /iK) 
comparable to that used in experiment 0- With these 
parameters, the equilibrium properties of the effective 
single-channel and BF model essentially coincide (for 
further details and values of parameters, we refer to 
Ref. 01). Therefore, to keep our discussion concise, 
we will focus on the single-channel theory noting that 
the parallel application to the BF model with appropri- 
ate physical parameters generates quantitatively similar 
results. At a field of ca. 1.0 gauss above the FR, the 
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FIG. 1: Variation of the chemical potential /j, with field B at 
T = for a mixture of fermionic 40 K atoms prepared in the 
(/ = 9/2, mj = -9/2) and (/ = 9/2, m/ = -7/2) Zeeman 
states at a density of n — 1.5 x 10 13 cm -3 (i.e. Tf = 0.35 (iK). 
The FR takes place at a field Bo = 202.1 gauss while the BEC- 
BCS crossover ([i = 0) occurs when B — Bo ~ 0.3 gauss. The 
inset indicates the scattering length a(B) = ab g (l — Jt% ) m 
the universal regime, where ab g is the background potential 
scattering length and AB the width of the resonance. Note 
that, in the present theory, we use finite range potentials with 
parameters which describe the FR also far outside the univer- 
sal region [ll|. 

condensate has an essentially BCS-like character while 
the experiment using the fast sweep technique observed 
the condensate starting from 0.5 gauss above Bq. To ex- 
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plore the entire region of interest, we choose as initial 
conditions field values B\ which span the entire crossover 
region (marked by stars in Fig. Starting from the 
ground state T = distribution, we follow the dynamics 
of the condensate after an abrupt switch to some different 
value of magnetic field Bp. 
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FIG. 2: Time dependence of the order parameter ampli- 
tude |Ak=o|/A oq following an abrupt switch from a field of 
(Si - B )/gauss = -0.5 (top), 0.0, 0.5, 1.0 (bottom) to 
(Bf — £>o)/gauss = —10, and from (Bi — Bo)/gauss = —0.5 
(top inset) and (Bi — Bo) /gauss = 0.5 (bottom inset) to 
(Bf — Bo)/gauss = —1.0. In addition to the constant phase 
velocity of Ak=o (found numerically to be set by the final 
state equilibrium chemical potential fj,), there is an additional 
time-dependent phase modulation whose characteristics mir- 
ror closely that of the amplitude oscillations. 

Figure [5] shows the time-evolution of the order param- 
eter |Ak=o|j normalised by the value, A eq (i.e. the value 
that it would acquire were the system to reach the T = 
ground state at the final field Bp). Here we have chosen 
a field Bp — Bq = —10 gauss deep within the BEC phase 
where the large binding energy of molecules allows their 
momentum distribution to be inferred from time of flight 
measurements 0- In contrast to the predictions of the 
constrained dynamics [||,IS0> these results show that (a) 
the coherent oscillations are substantially damped even 
at the level of mean-field, (b) the amplitude of the os- 
cillations is small, and (c) the order parameter |Ak=o| 



asymptotes to a value much less than the expected fi- 
nal state equilibrium value A eq . Referring to the insets 
of Fig. |3 one may note that, when the initial and final 
conditions are drawn closer, the period of the oscillations 
becomes longer and the effects of the damping more pro- 
nounced. Moreover, although the period of the oscilla- 
tions increases monotonically with A oq , the dependence 
is nonlinear and, referring to the bottom inset of Fig. [31 
one may note that A eq does not provide a ceiling for the 
magnitude of the oscillation. 

To interpret the generic behaviour, it is instructive to 
access the time-dependence of the pair wavefunction «k 
and the distribution function $k- Figure |3| shows that, 
when the abrupt switch takes place from B\ — Bq = 
0.5 gauss to Bp — Bq = —10 gauss, although there is a 
slight tendency to shift towards the final state equilib- 
rium distribution, the pair wavefunction and the density 
distribution remain essentially frozen close to the initial 
BCS-like distribution, exhibiting only small oscillations 
in time. In the absence of energy relaxational processes, 
the system is unable to significantly redistribute weight. 
By contrast, when the abrupt switch takes place from 
B\ — Bq = 0.5 gauss to Bp — B = —1.0 gauss, the prox- 
imity of the two phases allows the system to converge to a 
(non-stationary) modulated distribution whose (station- 
ary) envelope reflects more closely the final state equi- 
librium distribution. Referring to insets in the bottom 
panel one may note that a weakly perturbed conden- 
sate on the BCS side attenuates to its equilibrium value 
with the pair distribution showing only small modula- 
tions around the equilibrium in accord with the linear 
stability analysis of the weakly perturbed BCS system @ 
discussed in Ref. |l5j . In particular, one may note that 
here (and, indeed, for other values of the initial and final 
conditions) the hole-burning phenomenon predicted by 
the constrained dynamics does not appear. 

To assess the potential to observe coherent atom- 
molecule oscillations, one can monitor the time evolution 
of the number of condensed molecules, 



n mc (t) 



d 3 k K k (t) <fe(k,B F ) 



(3) 



Here (f>s,{k,Bp) denotes the wavefunction of the highest 
vibrational bound state being an exact eigenstate of the 
two-body problem Referring to Fig.0J one may note 
that the amplitude of oscillations is negligible (less than 
one percent). Although one may adjust the final field Bp 
to lie closer to the FR, the amplitude of the oscillations 
increases only slightly while the damping rate is also en- 
hanced. Moreover, the inclusion of processes beyond the 
mean-field approximation considered here would simply 
increase the damping rate and not enhance the oscilla- 
tions. One may therefore conclude that the observation 
of atom-molecule oscillations, following an abrupt change 
in the interaction strength, is infeasible for the entrance 
channel dominated resonances currently discussed in the 
context of the 40 K and 6 Li system. 
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FIG. 3: The pair wave-function \Kk\ (upper panel) and the 
density distribution function $k (lower panel) shown as func- 
tion of k = |k| with Bi — Bo = 0.5 gauss (dashed lines) 
and Bf — Bo = — 10 gauss after 50 fis (dashed-dotted line) 
and Bf — Bo = —1.0 gauss after 800 /is (solid line) following 
the abrupt switch as in Fig. The dotted lines signify the 
ground state equilibrium distributions at Bf — Bo — —1.0 
gauss included for comparison. The upper inset shows oscil- 
lations of |«k=o] for Bf — Bo = — 10 gauss. The lower two 
insets refer to a rather weak perturbation on the BCS side 
from Bi — Bo = 0.3 gauss (dashed line) to Bf — Bo = 0.2 
gauss. The dotted lines shows the equilibrium distribution 
while the solid line provides a snap-shot of the distribution at 
3 ms after the switch. Note that the harmonic modulations 
visible in the distribution functions translate to a single en- 
ergy scale of the same order of magnitude as the period of 
oscillations seen in |Ak=o|- 

To conclude, we have presented a numerical analysis of 
the dynamical mean-field equations for the single-channel 
theory of the FR following an abrupt field change. In 
the range of physical parameters appropriate to the 40 K 
system, the consideration of the two-channel BF (Bose- 
Fermi) theory does not change the results gwaniitatively. 
When applied to a theoretical regime where the pop- 



ulation of the closed channel states below resonance is 
high, the numerical findings do not change qualitatively. 
Relying on the deterministic time-evolution of the ini- 
tial stat e according to the Hciscnberg mean-field equ a- 
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FIG. 4: Time-dependence of the relative number of condensed 
molecules (n mc (t) — n mc (0))/n mc (0). Here we have used the 
same field values as that used in Fig. |21 (main) with (Bi — 
Bo)/gauss = —0.5 (top), 0.0, 0.5, and 1.0 (bottom) and (Bf — 
Bo)/gauss = —10. For these field values, when normalised to 
one half of the total atomic density, n mc (0) = 0.0005, 0.012, 
0.1, and 0.3 respectively. 



tions of motion, our results differ substantially from the 
findings of the constrained dynamics considered previ- 
ously [1, IE • Over a wide range of initial conditions, 
we observe substantially damped oscillations with an am- 
plitude strongly dependent on initial conditions and a fre- 
quency set by the interaction strength after the switch. 
To assess the capacity for BCS-BEC like atom-molecule 
oscillations following an abrupt change in the interac- 
tion, we have chosen initial and final conditions to span 
the crossover from the BCS to the BEC limits. We have 
found that the amplitude of atom-molecule oscillations 
is negligible and the distribution is essentially frozen to 
the initial, a behaviour reminiscent of an orthogonality 
catastrophe. We conclude that, in the entrance channel 
dominated resonances observed in 40 K and 6 Li, the ob- 
servation of atom-molecule oscillations is infeasible. The 
rigidity of the condensate wavefunction distribution at 
short-time scales, where the time evolution is mean-field 
in character, supports the method of the fast sweep as a 
reliable technique to probe the fermionic pair condensate. 
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